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Ch ■ Abstract 

We study the existence and properties of wormhole throats in modified f{R) 
gravity theory. Specifically, we concentrate on the cases where the lapse is not 
! necessarily constant, and hence are not limited to the zero tidal force scenarios. 

In the class of theories whose actions are generated by Lagrangians of the form 
Q I f{R) = J2'^n we find parameters which allow for the existence of energy 

condition respecting throats, which do not exist in Einstein gravity. We also 
^ \ consider the effect of the modified action on the anisotropy of the models, and find 

bX) that modified gravity can minimize the amount of anisotropy required to support 

the existence of a throat. In both these respects, the sector containing theories 
^SJ ■ with positive n is more promising than the negative n sector in comparison to 

^ \ Einstein gravity alone, with large n being most favorable. 

o 

. PACS numbers: 04.50.Kd 04.20.Gz 
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1 INTRODUCTION 

General relativity is arguably the most successful theory of gravity to date. The theory, 
^ \ aside from being pleasing on grounds such as diffeomorphism invariance, has passed 

all major tests in the solar system. On scales larger than the solar system, general 
relativity, supplemented with non-standard matter fields such as dark matter or dark 
energy, explains the dynamics of the large-scale structure of the universe rather well. 

However, there has long been an interest in extended theories of gravitation, which 
possess general relativity in some limit (for example, see [1], [2].) These extensions 
are mainly motivated by the desire to eliminate the potentially "exotic" non-standard 
matter fields mentioned [3], [1], or to mimic some low energy quantum gravity effects 
which are thought to manifest at very high curvatures [5], [6]. Other motivations 
are perhaps more academic; we have only been able to study gravity effectively in 
its weak limit, and hence cannot be certain that Einstein gravity holds in stronger 
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gravitational fields. Therefore it is of interest to study gravitational theories which are 
diffeomorphism invariant and give Einstein gravity in an appropriate limit, but deviate 
from Einstein gravity in some way outside of the realm where gravitational effects have 
commonly been observed. This, admittedly, gives us much freedom regarding the types 
of theories that may be successful candidates. 

Of particular interest are the class of theories known as f{R) gravity. These 
theories consider gravitational fields generated by actions of the form [7j: 

S=^j ^f{R)d^x + j ^C,Ug)d^x, (1) 

where f{R) is some function of the Ricci scalar, R and Cn^^t{g) is the matter Lagrangian 
density. The equations of motion generated via varying this action with respect to 
the metric, and demanding that this variation vanish. These equations of motion are: 

WR'^ - IfiR) - V^V,^ + ^/^^ = , (2) 

where as usual is the stress-energy tensor of the matter field which comes from 
varying the matter action. 

A number of earlier generalizations to Einstein gravity fall within the paradigm of 
f{R) theories. For example, theories with f{R) = R + a2R^ have long been studied 
(called "/^-squared" gravity. See, for example, [8]-[Il]). The Starobinsky inflationary 
theory p!5] may be the most popular application of i?-squared gravity theory. In the 
opposite regime, where curvature is low, there has been less interest in modifying 
Einstein gravity. However, it is worth mentioning that f{R) theory with inverse 
powers of the Ricci scalar has been considered as a possible candidate to explain the 
observed accelerated expansion of the universe |16]-|24j. Inverse Ricci terms may also 
appear in certain sectors of string/M theory (see |25) . |26j for example). The vacuum 
state for inverse- i? theories is not Minkowski space-time, but is instead either deSitter 
or anti-deSitter space-time. Admittedly, these inverse theories are highly constrained 
from stability considerations or solar system tests but may not be completely ruled 
out. (See ^27j-|36j and references therein for a discussion of the restrictions.) For 
generality we also include this sector in our study. A thorough study of solutions in 
i?" cosmology may be found in [37j as well as [38] and a complete discussion of f{R) 
modifications to aid acceleration may be found in [3^. Outside of cosmology, and 
perhaps more relevant here, studies have been performed that have more bearing to 
stellar physics |40j-|48j. In another popular extension, namely that of f{T) gravity 
where T is the torsion scalar, such solutions have also been studied [l9], [50] . 

Given that integer powers of R have arguably generated the most interest in f{R) 
extensions, we propose here to study the properties of wormhole throats in gravity 
theories given by 

f{R) = Y,anR'', (3) 

n 

with the an constants. This is quite general as any f{R) analytic in R may be 
expanded as such a power series (in the positive n sector) such as exponential gravity 
[51] . For completion, we allow for both positive and negative powers of R in the 
action (hence extending to f{R) expandable in negative powers as well). Of course, 
this series includes the n = 1 term (Einstein term) and can include n = (cosmological 
constant) . 
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The wormhole is of interest as it may provide a simple model for the space-time 
foam thought to be manifest at very high energies, where the corrections to Einstein 
gravity due to higher curvature terms may be important in encompassing some quan- 
tum gravity effects. As well, exhaustive studies have been performed on wormhole 
geometries in Einstein gravity and their properties (mainly in the static cases with 
spherical and axial symmetry) are now well known in Einstein gravity. For example, 
it is an interesting property that static wormhole throats necessarily must violate the 
weak/null energy conditions in Einstein gravity [52], [53] and, furthermore, must be 
anisotropic. It has been of much interest in wormhole physics to find either some- 
what realistic matter models which can meet these slightly exotic conditions in a way 
required to support a wormhole throat, or else to consider alternative gravitational 
theories which may allow for energy condition respecting matter. In the latter vein, 
due to the complexity of the resulting equations, the zero tidal force class of wormholes 
is most often studied {gu = const.). Even then, analytic models are hard to come by 
and one often resorts to numerics. The zero tidal force models are interesting due to 
their tractability, which allows one to study important properties of these geometries. 
However, they may not be very realistic from a physics perspective though, as they 
yield a constant frequency shift in inhomogeneous structures possessing a preferred 
center and, in the weak-field limit, a constant Newtonian potential. 

In this work we study wormhole throats and without the restriction of zero tidal 
force, although we do consider the zero tidal force cases as well for completion. The 
consideration of the throat is in many ways more general than considering the asymp- 
totics, as the most salient features of a wormhole occur in the throat region. For 
example, in Einstein gravity, the necessary violation of energy conditions in static 
spherically-symmetric wormholes occurs in the neighborhood of the throat, regardless 
of asymptotics. In fact, in the Einstein gravity scenarios, if the matter field falls off 
sufficiently fast, or is patched to a vacuum or other solution, the properties far from 
the throat do not generally differ greatly from similar systems with trivial topology. 
As well, it has been argued (for example in [53]) that the global topology is too limited 
a tool to study wormholes, and a local geometric analysis near the throat is generally 
more useful in discerning interesting properties of wormholes. In such cases, the throat 
does not necessarily coincide with the common definition of a wormhole, which relies 
on global properties, and is viewed as an interesting object in its own right, capable of 
describing more general scenarios than just wormholes. We therefore now concentrate 
on the near throat region and study the properties in this vicinity. Admittedly, in 
some cases, demanding flatness at infinity would place extra restrictions on the prop- 
erties of the wormhole [55j ) [55] but as the asymptotics of the universe are not exactly 
known, and we are interested in the existence of throats only, we will not consider 
such restrictions here. 

In section 2 we briefly review some background material and present a non- 
standard coordinate gauge that is more suited for wormhole analysis than the standard 
spherical coordinate chart. We discuss the mathematical construction of the throat in 
this chart. This is followed by a study, which is analytic when possible but otherwise 
numerical, of various scenarios depending on the terms present in the gravitational 
action. Of special importance, due to their physical relevance, are the non-zero tidal 
force models, which we examine in some detail. We particularly concentrate on the 
properties of energy conditions as well as the degree of anisotropy. In general, we 
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find that the parameter space of Einstein gravity supplemented with terms where n 
is positive is more favorable in both these respects than Einstein gravity alone for 
a large range of parameters. Supplements with negative n have either a very small 
parameter space where they improve the energy conditions or else tend to worsen the 
situation when compared to Einstein gravity alone. Finally, in section 3, we conclude 
the study. 



2 THE MODELS 



As is common in wormhole studies in Einstein gravity, we will utilize here an anisotropic 
fluid source whose stress-energy tensor is given by: 

T^u = {P+ Pt)u^Uu + Pt K + iPr - Pt)s''s^ . (4) 

Here p, pt and Pr are the energy density, the perpendicular (to the inhomogeneous 
direction) pressure, and the parallel pressure respectively as measured in the fluid 
element's rest frame. The vector u'^ is the fluid 4-velocity and s^^ is a space-like vector 
orthogonal to u^. These vectors satisfy: 



% = -l, s% = +l, ^^% = 0. (5) 



Since we are interested in the properties of the actual material generating the gravi- 
tational field, we do not transform the system into an effective scalar-tensor system 
but instead keep the geometry and the material properties explicit. 
We impose spherical symmetry and hence may write the metric as 

ds^ = - e*('-) + e-^'-) dr^ + dO^ + sin^ (9) dip" , (6) 

= - dt' + [l + [drP{r)f] dr^ + dO'' + W{e) V , 

where the function P{r) describes the profile of the three geometry for fixed 9 and (p 
(see below for details). In this chart the Ricci scalar is given by 



R 



2r2 



1 + {P') 



2 2 



A{P'f - 4r(P')^*' - 4r$' + SrP'P" + 2r^P'P"^' 



+ 4{P'f - 2r\P')H" - r\P')\^'f - 2r^^" - r\^'f , (7) 



where the prime denotes differentiation with respect to r. Although the above is 
the most common form of metric system to locally describe spherically symmetric 
gravitational fields, this coordinate system is not ideal for studies of wormhole throats 
and hence we use a different chart, which we now present. 



2.1 WORMHOLES 

We consider here static throats, which must be anisotropic in Einstein gravity, in the 
class of modified gravities given by f(R) = ^OnR"' for both positive and negative 

n 

n. As mentioned previously, it is common in the literature on wormholes to consider 
zero tidal force models via the condition guir) = constant. However here we relax 
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this restriction and study more general scenarios. As well, although the spherical 
coordinate chart of ([6]) can be utilized for wormhole throats if care is taken, it is 
not optimal, as grr{r) — )■ oo as r — )■ (see Figure [J and equation and the 

most interesting properties of wormholes are arguably found near the throat region. 
Therefore, we use a different chart for the study of wormholes which essentially involves 
tilting the standard spherical chart by 7r/2. The wormhole throat is then given via the 
creation of a surface of revolution of the profile curve, r = Q{x), as shown in Figured] 
(also see figure caption). In this new chart, the space-time metric's line element may 
be written as 

ds2 = .gi-W ^^2^1^^ [d^Q{x)f^ dx^ + Q\x) de^ + Q\x) sin^ie) V • (8) 

Note that in this new chart the metric is analytic at the throat since dxQ{x) — )• 0, 
and hence gxx{x) — )• 1 as one approaches the throat (x = o|3- Also, only one chart 
is required now to cover the wormhole, as opposed to two charts in the standard 
spherical coordinates. The function Q{x) must possess the following properties: 

i) Qo := Q(0) > 0, 

ii) Q[, := Q'{x)\x=o = 0, 

iii) Q"{x) > in some neighborhood of the throatH. 

Aside from the above properties we make the mild assumption that Q{x) is analytic 
in some nonzero domain about x = 0. 

In the above coordinate chart, the Ricci scalar is given by 



R 



1 



2Q2 1 + (Q' 



,2 2 



4 - 4Q(Q')^' + 4(Q') - 8QQ" - 4QQ'«>' + 2Q'Q'Q"i^' 



where here and subsequently the prime denotes differentiation with respect to x. We 
do not expect any serious pathologies (such as infinite tidal forces) in our analysis as, 
in the domain of validity of the coordinate chart (which is larger than the "r" chart, 
and covers the throat) the metric and its derivatives are well behaved and hence we 
expect the orthonormal Riemann tensor components to also be well behaved. We 
show later that the function Q{x) describing the spatial geometry, is C°° and that 
^>(x) is also very smooth. 

One issue which is of great interest in wormhole physics is the study of energy 
condition violation of the matter supporting the wormhole. It is known that in Einstein 
gravity energy conditions must be violated by static wormholes somewhere in the 
vicinity of the throat |52], [53] . It is possible that in more complicated gravitational 
theories, one may circumvent this issue and possess a throat region which respects 
energy conditions. An analysis in this vein was performed in |56) where the analysis 



^As an aside, in this coordinate system the exterior Schwarzschild metric is given by Q{x) — 

2 

2M + -§7-!. The horizon is located at a; = 0. 

More precisely, if Q's first non-zero derivative (higher than first order) at a; = is of even order, 
the function attains a local minimum if this derivative is positive, and hence we have a wormhole 
throat. If its first non-zero derivative is of odd order, it is a point of infiection and therefore does not 
describe a wormhole throat. 
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Figure 1: a) In the standard spherical chart the wormhole throat is generated by the 
surface of revolution created by rotating the profile curve x = P{r) — P~{r) IJ P+{r) but 
this leads to a metric singularity at the throat (r = ro) where drP(r) — oo. A new chart 
is created in b) via rotating the coordinates by tt/2, and the surface of revolution (inset) 
is generated by rotating the curve r = Q{x) = P~^{x). This does not have the metric 
singularity at the throat (a; = 0). Also, a single chart can now cover both regions. 



was limited to n = — 1 and n = 2 extensions only. A very complete analysis may 
be found in [57], [58]. In the above cited studies, however, due to the complications 
involved, the zero tidal force assumption was made. This simplification eliminates 
some of the complications but still allows one to glean a number of the interesting 
properties that f{R) wormholes possess. However, constant redshift, and constant 
weak-field Newtonian potential, from various regions of an inhomogeneous gravitating 
object may not be particularly realistic and we wish to remove this assumption and 
allow for a non constant ^(x) function. This can alter the physical properties of 
the throat significantly, as we shall see below. Aside from yielding more complicated 
equations, a non-constant <I>(x) now presents us with the problem of how to prescribe 
this function. For this we appeal to some physical considerations. It is more physical 
to demand some realistic properties on the matter fields than it is to blindly prescribe 
although the former is much more difficult than the latter, as one needs to solve 
a complicated differential equation for <I>(x). With this in consideration, what we 
initially do is prescribe some reasonable energy density profiles, p{x), for the matter 
field (eg. large positive energy density with zero slope at the throat and monotonically 
decreasing outward towards the surface or infinity, using parameters where this is 
possible) and numerically solve for $(x) using the numerical code COLSYS [59] . This 
will then give us an idea of what a realistic ^(x) function should look like, and we 
use similar functions for subsequent studies. (The spatial geometry must also be 
prescribed, and this poses no problem as we apriori know that, spatially, we are dealing 
with a spherically symmetric wormhole.) We simply use this method to provide an 
idea of what the functions <I>(x) should look like in realistic scenarios and then use 
similar functions for our studies. It should be noted that the function <I>(x) controls 
the presence of event horizons. We avoid horizons by ensuring that ^(x) does not 
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approach — oo anywhere in the domain of study. 

In the following analysis wc wish to discern how the various terms in the action 
affect the wormhole solution. In the case of f{R) = '^anBP'-, each term in the sum 
contributes independently to the left-hand side of the field equations, and hence also 
contribute separately in determining the properties of T^. Therefore, it is sufficient 
to study individual terms in the power series of f{R) separately to see how they each 
contribute individually to the matter field. Note that in this approach, a„ is simply 
an overall multiplicative constant, and hence we set it equal to one. One may see 
the effect of changing a„ by simply rescaling the results below by whatever value of 
an one chooses (with the exception of the anisotropy studies, for which we study the 
effects of varying a„). Negative values of a„ also reflect the graphs about the z = 
plane. 



2.1.1 n=l 



This case, of course, corresponds to Einstein gravity. It is worthwhile presenting this to 
summarize the results of Einstein gravity in this coordinate chart, and for comparison 
with other powers of n later. Regardless of the powers present in the full action, it is 
expected to have an n = 1 term. We summarize the n = 1 results as follows: 



p :=Kp = —kT^q 



1 + iQ'Y - 2QQ" 



1 + {Q'y 



2 ' 



1 



(lOi) 
(lOii) 



4g 



1 + {Q'Y 



-2QQ'Q"^' + 2Q{Q'fV + Q{^'f + 2Q$"] 



(lOiii) 



There are also the following combinations which occur in the energy conditions: 



P + Pr 



P+Pt 



Q[i+(Q'r 
1 



(Hi) 



4g2 



1 + iQ') 



[4 + 4{Q'Y - AQQ" + 2QQ%' + 2Q{Q'f^' 



^Q^(Cl'f{^'f - 2Q2q'q"$' + 2Q2(q')2$" + Q2($')2 + 2Q2$"] . (iiii) 



The above expressions, though not overly complicated, do not shed much insight 
into the behavior of the matter fields near the throat, where we are most interested 
in. We therefore expand the relevant expressions in Taylor series about the throat 
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(x = 0): 



1 O'" 1 

~P=ir[l- 2Q,Ql] - 2^x + [4Ql{Q'^f - Ql - ggQ^] + 0{x^) , (12i) 



P + Pr = -2^ + [Q[;$'o - 2Qo'] ^ + 4 4Qo(Qo)' - QoQo' + ^QoQo"<J>: 



1 



2 



'0 



+ QoQon + (Qo)'J + 0{x-') , (12ii) 

+ QoK - 2Q'o" + Q'^^'o] X + . (12iii) 

where the zero subscript indicates that the quantity is evaluated at x = 0. 

Prom the ()12ip . it can be seen that the energy density in the throat region may be 
made positive. Regarding ()12iip . if Q'o is non-zero it must be positive (as the throat 
is a local minimum) and (I12iip must therefore be negative near the throat. If Qq is 
zero, then the condition for a local minimum implies that Qq is also zero, and hence 
the lowest order term which would contribute near the throat is the fourth derivative 
term, Qq" ■ Since this fourth derivative must then be positive under the condition of a 
minimum, and it appears with a negative sign in p2iip . this contribution is negative. 
In such a scenario energy conditions are (barely) met at the throat, but are violated 
as one moves away from the throat. Similar arguments apply to higher derivatives in 
case the fourth derivative vanishes at the throat (although one needs to study higher 
order terms in the expansion, which for brevity we did not write). It is in this way 
that energy conditions must be violated in the vicinity of a static wormhole throat in 
Einstein gravity. 

For comparison with some of the modified gravity results to follow, which, due to 
the complicated expressions they yield, must be studied via computational methods, 
we choose a spatial geometry governed by 

Q{x) = Acosh^|-^ , (13) 

as this function possesses all the salient properties to describe a throat. The parameter 
A represents the radius of the throat and xq represents the degree of "flare-out" of 
the wormhole. Both these parameters are strongly related to the degree of energy 
condition violation, and hence we pay particular attention to these quantities. For all 
cases studied in this work, when an explicit form of Q{x) is required, we use the form 
in ([13]). 

For the function $(x) we appeal to physical considerations (with the aid of the 
numerical code COLSYS, as mentioned previously). For positive n we expect ^>(x) to 
slowly asymptote to a constant value far away from the throat, where the geometry 
is expected to approach Minkowski space-time. For negative n (and for generality, for 
all n) the asymptotic value should approach the deSitter or anti-deSitter value. As 
we are interested in the near-throat region only it is not crucial that the asymptote 
is manifest in the domain of consideration, but an indication of such an asymptote is 
desirable. We consider both scenarios where gu = — e*^^) smoothly increases towards 
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the asymptote, as well as scenarios where gu smoothly decreases towards the asymp- 
tote. (One cannot necessarily rule out one scenario over the other in extended gravity 
theories.) Specifically, we choose fitting functions of the form 



gtt{x) =Ao 



1 + x2 



Co 



for concave-down gu , 



(14a) 



9tt{x) 



Aq + Bqx^ 

Co + x2 



for concave-up gu , 



(14b) 



Where Aq, Bq and Co are fitting constants to fit generally to the COLSYS numerical 
results. Since the wormhole profile chosen in (jl3p is symmetric about the throat, we 
also make similar symmetry demands on gu- There is no great loss of generality in 
doing this since if one wishes to model a non-symmetric throat, the solution presented 
can be viewed as being valid only on one side of the throat, and a different solution 
can be patched to the other side. The forms of (jl4al b) are used throughout this 
manuscript in the analysis of non-zero tidal force models. 

In the figures below (Figures [2] -[7|) we display the results of Einstein theory (n = 1) 
for the spatial geometry governed by (jl3p . The first set of figures (figs. [2^-d) show 
the behavior of gu, along with the energy conditions (jlOip and (jllil ii) as a function of 
the flare-out parameter, xq. Although the energy density is positive (as is p + Pt), the 
other energy condition is not positive. We can also see here the well-known situation 
that the greater the degree of fiare-out (equivalent to small xq) the more severe the 
energy condition violation in Einstein gravity. 

Figure [3] displays the behavior of gu , along with the energy conditions (jlOip and 
(lllit ii) as a function of the throat radius, A. Note that although the energy density 
can be made positive near the throat if the throat is not too large, the quantity (jllip 
is negative in the throat vicinity, confirming the analytic analysis above. 

In Figure S] we display the anisotropy, pt — Pr, as a function of throat radius. Note 
that larger throats require less anisotropy (i.e. are "more isotropic"). 
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(c) p + pr (d) P + Pt 

Figure 2: Einstein gravity, gtt decreasing, throat radius=0.05, varying flare-out xo- 
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Figure 4: Anisotropy of throat region in Einstein case, gu decreasing, xq = 1, varying 
throat radius. 
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In Figures O [6] and [7] we present a similar analysis as above, except that in these 
figures, gu is concave-up. Note the similarity of these results to the previous analysis. 




Figure 5: Einstein gravity, gu increasing, tliroat radius=0.05, varying flare-out xq. 
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0.1 1.2 

tt(x) (b) p{x) 

+ Pr (d) p + 

1 gravity, increasing^ a;o = 1, varying throat radius. 


-0 

Figure 7: Anisotropy of i 
throat radius. 


0.1 1.2 

hroat region in Einstein case, gtt increasing^ a;o = 1, varying 



2.1.2 n = 2 

We now study the n = 2 contribution. Adding an term to the gravitational 

action is arguably the most popular supplement to the Einstein-Hilbert action (see 
[8j-|14j). For example, this modification has been utilized to drive inflation purely 
from the gravitational sector (the Starobinsky inflationary theory [15]) or to capture 
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some low-order quantum corrections [5]. 

As mentioned earlier, since in the theory given by the action ([T|) each term in the 
series contributes to the energy conditions separately, it is more instructive to look 
at this contribution on its own, to see whether it helps or hinders with respect to the 
energy conditions. For the anisotropy analysis, we add the Einstein term to it as well 
to see if the modification can make the system more or less isotropic compared to 
Einstein gravity alone. We consider zero and non-zero tidal force cases. 

Zero tidal force: In the case of zero tidal force (<I>(x) = const.) we can get a handle 
on the behavior of the energy conditions at the throat by performing a series expansion 
as was done in the Einstein case. The results are summarized as: 

Wo 

+ ^ [QoQo" - 25Qom'Qo + ^QoQo] X + 0{x^) , (15i) 

P + Pr = ^ [QoQo' - ^QoiQ'of + 2iQ'^?] 

+ ^ [QoQo" - 25Qo(Q[)')'Qo" + X + 0{x^) , (15ii) 

~P + Pt = ^[^ + 2Qg(Qo)' - 3QoQ[)'] + ^Q'o" (2goQo - 1) ^ + • (ISiii) 

For 02 > 0, which is the more physical sector [1], it can be seen that the above can 
all be made positive at the throat (for example, by considering models where Qq = 
and Qq' > 0). By analyticity in a non-zero neighborhood, in such scenarios there 
therefore must be a non-zero domain about the throat for which the functions are 
non-negative. Therefore, in this class of models, throats (and we stress again here 
that we are not considering asymptotics at infinity) are allowed which respect energy 
conditions. This will remain true even when adding the Einstein term to the action, 
as in the Einstein case the near-throat violation of energy conditions may be made 
arbitrarily small [60], [61] independently of the a2 parameter. 

Non-zero tidal force: Unfortunately, for non-zero tidal force the analytic expressions 
with a general Q{x) and <I>(x) are very long and complicated, even as a near throat 
expansion, and not very revealing. We must therefore specify these functions in order 
to perform numerical studies. For this purpose we choose the same function as in 
the Einstein case so that comparisons may be easily made. That is, we choose the 
profile Q{x) as given in (fT3|) since, to reiterate, this function possesses all the required 
properties to describe a throat. Recall that the parameter A represents the radius of 
the throat and xq represents the degree of "flare-out" of the wormhole. We also use 
similar <I>(x) functions as presented in the Einstein case. 

As mentioned previously, we set Q2 = 1, as all results here can be rescaled by 
whatever value of 02 one wishes to study, including negative values, which result in a 
reflection about the horizontal planes in the graphs. As mentioned above, it should 
be noted that a2 > is the preferred model. The flrst set of results are summarized 
in Figure [8] (also please refer to figure captions for details). 
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Figure 8: contribution, gu decreasing, throat radius=0.05, and varying a;o. 



Note that like the zero tidal-force case, there is a large region of the parameter 
space where all energy conditions can be respected. In fact, for positive 02, a large 
flare-out (small xq) actually helps the condition p + Pr > near the throat. 

Also of interest is the anisotropy. Namely, can the presence of the extra terms in 
f{R) gravity lessen the amount of anisotropy required to support a throat, compared 
to Einstein gravity alone. To study this we present Figure [9] where the anisotropy, 
defined here as Pt—Pr-, is presented for the Lagrangian i? + a2-R^ and is plotted in the 
vicinity of the minimum anisotropy curve. Note from this figure that anisotropy is a 
minimum for a2 < 0, and not for the Einstein case (02 = 0), although the minimum 
curve is very close to 02 = 0. 

In Figure [To] we present another non-zero tidal force analysis (with the same gtt{x) 
as in the previous case) where, instead of allowing xq to vary, we vary the throat radius. 
This study is interesting as even in Einstein gravity the size of the throat affects the 
amount of energy condition violation. We find that, generally, smaller throat radius is 
more favorable for respecting energy conditions. (This result is reversed for 02 < 0.) 

Next we consider non-zero tidal force models where gtt{x) is increasing instead 
of decreasing. The figures below (figs. [H] and [T2|) summarize the results for varying 
Xq and 02. Note that in Figure [11] the graphs mimic the previous scenario with 
decreasing guix), indicating insensitivity to the form of the lapse function. Therefore, 
in this scenario, again energy conditions can be respected for the more physical sector 
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Figure 9: Minimum anisotropy region: R + a2R^ contribution, gtt decreasing, throat 
radius=0.05, xq ~ 1, and varying a2- 




(a) p{x) (b) p + pr (c) p + Pt 



Figure 10: contribution for xq — 1, gtt decreasing and varying throat radius. For 
negative a2 the graphs should be flipped around the plane z = 0. 



of positive 02- We also present the anisotropy in Figure [12] for the full Lagrangian of 
R + a2-R^. Anisotropy is again minimized for Q2 7^ and in the negative sector. 

Finally for this sub-section, we present the results for the increasing gtt{x) but 
where the throat radius is allowed to vary. These results are summarized in Figure [T3j 
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Figure 11: contribution, g^^ increasing, tliroat radius=0.05, and varying a;o. 




0.06 -0.0002 



Figure 12: Minimum anisotropy region: R + a2R? contribution, gtt increasing, tliroat 
radius=0.05, xq = 1, and varying a2- 
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(a) p{x) 



(b) p + pr 



(c) p + vt 



Figure 13: i?^ contribution for = 1, gtt increasing, and varying throat radius. For 
negative a2 the graphs should be flipped around the plane z = 0. 



Note that again smaller throat radius is more favorable for energy conditions than 
larger throat radius if 02 is positive. (The result is reversed for negative 02.) Param- 
eters exist where energy conditions may again be satisfied, even when adding these 
results to the energy condition violating Einstein gravity terms. 

2.1.3 n = 3 

The n = 3 scenarios are also of interest and it has been shown in a cosmological 
setting how this case can be transformed to a system consisting of Einstein gravity 
and a scalar field with a self-coupling proportional to ast/)^ [62j. 

Zero tidal force: Again for the ^{x) = const, scenario we present an analytic expan- 
sion about the throat, although, due to the complexity of the coefficients, we produce 
only the lowest order term. The results are: 

p = ^ [1 - 24Q^g^'Q'o'" + 6QoQ'o' - 24g^(Q^")' + I^QIQ'^;' - 12g2(Q^')2 
^0 

-^QQliQlf + 96g^(Q^')1 + 0{x) , (16i) 

~P + Pr = '^-^ [Q'o - ^QlQoQ'o - ^QliQo? + -^QlQo - i2Qg(Q'o')' 

+lQQMf] + 0{x) , (16ii) 

~P + Pt=^^[i- 5QoQ;,' + ^QliQ'l? - ^QliQlf] + 0{x) . (I6iii) 

Here, in the case of positive as, ()16ip - (|16iiip can in principle be made positive. 
The simplest models obeying energy conditions would again be those where Qq = 
(which, due to the condition for a minimum, also requires Qq = 0) and Qq" > 0, 
as was the case for the corresponding situation in n = 2. Note that if Qq = 0, the 
contributions from Einstein gravity, from ()12ip - ()12iiip . exactly at the throat are 
non-negative (although ()12iip becomes negative in some neighborhood away from the 
throat), and hence with the addition of these a^B? contributions, f{R) = R + a'^B? 
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gravity can be made to obey energy conditions (and therefore, including the previous 
analysis, f{R) = R+a2R'^+a3R^ gravity can also be made to obey energy conditions). 



Non-zero tidal force: The analysis of the non-zero tidal forces proceeds in a similar 
order as the n = 2 case. Due to the interest in keeping the paper of reasonable length, 
we present all the graphs (figs. [U]- [T9|) for this case first, and then briefly comment 
on the results afterward in table 1. (Please refer to figure captions for details.) 




(a) gtt{x) 





(c) p + Pr 



W p + pt 




Figure 15: Minimum anisotropy region: R + a^R^ contribution, gtt decreasing, throat 
radius=0.05, xq = 1, and varying a^. 
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(a) p{x) 




(b) p + pr 




(c) P + Pt 



Figure 16: contribution for xq = 1, gtt decreasing, and varying throat radius. For 
negative the graphs should be flipped around the plane z = 0. 




Figure 17: contribution, gu increasing, throat radius=0.05, and varying xq. 
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Figure 18: Minimum anisotropy region: R + a^B? contribution, increasing, throat 
radius=0.05, xq = 1, and varying as. 




Figure 19: contribution for xo = 1, gu increasing, and varying tliroat radius. For 
negative as the graphs should be flipped around the plane z = 0. 
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We summarize the results for n = 3 in table 1. 



Table 1: Summary of a^iB? contributions. 



Parameter 
studied 


zero tidal-force 


gtt{x) concave-down 


gtt{x) concave-up 




Possible to make energy 
conditions positive. 


Tendency to make pos- 
itive contribution to all 
energy conditions for all 
range of studied near 
the throat for positive 

r\i n fi n H n Pccfi T\rp p(~jiTt"T"i — 

bution for negative a^. 


Tendency to make posi- 
tive contribution to all 
energy conditions for 
positive as, and negative 
contribution for negative 

L-t3 ■ kJiiU W o iiioV^iioi Lri V i U V 

to form of gtt {x) . 


Throat radius 


Possible to make energy 
conditions positive. 


Energy conditions more 
positive for small throat 
radius (with as > 0. Re- 
verse reported results for 
as < 0). 


Energy conditions more 
positive for small throat 
radius (with as > 0. Re- 
verse reported results for 
as < 0). 


Minimum 
anisotropy for 
R + a^R^ 


Not studied. 


q;3 < minimizes 
anisotropy. 


as < minimizes 
anisotropy. 



2.1.4 n = -l 

Zero tidal force: Again we begin with an analytic analysis of the zero tidal-force 
model. The results are summarized as: 



QoQo - tj 



1"^2 

4(2Qc 

-24Qg(Q^')3 _ 8Q^g'-Q^' + 32Q^(Q'^')4] + o{x) , (171) 



P + Pr= ,"-'T' ,,4 [3Qo + ^QW + l2QUQ'^f - ^QM? 
I {2QoQq - i) 

-AQM'' - ^QlQ'^Q'^' + iQQliQo)'] + 0{x) , (I7ii) 
P + Pt = ,.,^~'T' ,,2 [QoQo - 1] + 0{x) . (17111) 

Here it can be seen that with inverse powers of i? a peculiar singularity occurs. The 
singularity occurs when 2QQ" = 1 and corresponds to the vanishing of the Ricci scalar 
at these points. This is not unexpected and it is not a curvature singularity, but does 
herald a problem with the equations of motion when R = 0. One common way to 
attempt remedy this situation in studies of inverse- i? gravity is to postulate that the 
(gravitational) vacuum state of the theory is not Minkowski space-time, but is instead 
deSitter or anti-deSitter space-time |20j . However, when one is not far away from 
sources (as in the study here) there can be curves or surfaces in the space-time on 



22 



which R = and these must be excluded. Hence the parameter space studied here 
will not include parameters near these pathologies, which we excise. 

Here the vanishing of Qq does not yield an energy condition respecting throat, 
unlike in the positive n cases. (Either p + pt < for a_i > or else p + Pr and p < 
for a_i < 0.) Note that in principle it may be possible to respect energy conditions 
near the throat. One can say that if QqQq > 1 and furthermore if Qq" is large then it 
may be possible to respect energy conditions, but the situation is not obvious. Hence 
here the numerical results are needed, which we present for the more general non-zero 
tidal force scenarios. 

Non-zero tidal force: If one is considering the region far from the throat, then gtt{x) 
should, as mentioned above, asymptote to deSitter or anti-deSitter space-time. Our 
choices for gtt{x) can accommodate this, but we are only interested in the near-throat 
region, so strictly speaking it is not a requirement. 

Again, due to the interest in keeping the length reasonable, we present all the 
graphs for this case first (figs. [20] - [251) . and then briefly comment on the results 
afterward in a summary table (table 2). 




(c) p + Pr (d) P + Pt 



Figure 20: R ^ contribution, gtt decreasing, throat radius=0.05, and varying xq. 
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Figure 21: Anisotropy of throat region: R + a-iR ^ contribution, gtt decreasing, throat 
radius=0.05, a;o = 1, and varying 




Figure 22: contribution for a;o = 1, gtt decreasing, and varying throat radius. For 

negative a-i the graphs should be flipped around the plane z = 0. The region around 
the singularity has been omitted. 
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Figure 23: R ^ contribution, increasing, tliroat radius=0.05, and varying xq. 




Figure 24: Anisotropy of tliroat region: R + a-iR ^ contribution, gu increasing, tliroat 
radius=0.05, xq = 1, and varying a_i. 
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(a) p{x) (b) p + pr (c) p + vt 



Figure 25: R^^ contribution for = 1, gu increasing, and varying throat radius. For 
negative Q!_i the graphs should be flipped around the plane z = 0. The region around 
the singularity has been omitted. 



We summarize the results for n = — 1 in table 2. 



Table 2: Summary of a_ii? ^ contributions. 



Parameter 
studied 


zero tidal-force 


gtt{x) concave-down 


gtt{x) concave-up 




Respecting energy condi- 
tions seems difficult at 
best. 


Little sensitivity very 
near the throat on 
variations of xn- 


Little sensitivity very 
near the throat on 
variations of xo . 


Throat radius 


Respecting energy condi- 
tions seems difficult at 
best. 


Possible to respect en- 
ergy conditions near sin- 
gular region for negative 
a_i for a small section of 
parameter space. (Only 
on one side of the singu- 
lar region.) 


Seems possible to make 
small positive contribu- 
tion to all energy condi- 
tions (see general n sec- 
tion) for very small re- 
gion of parameter space. 


Minimum 
anisotropy for 


Not studied. 


Magnitude of anisotropy 
larger than correspond- 
ing pure Einstein case re- 
gardless of q;_i. 


Magnitude of anisotropy 
larger than correspond- 
ing pure Einstein case re- 
gardless of q;_i. 



2.1.5 n = -2 

Finally we present here the contribution from a_2^~^- Much of this scenario mimics 
the a-iR~^ contribution and hence we also summarize the results in a table (table 3) 
after all the graphs (figs.[26l- [3T]) . 
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Zero tidal force: At the throat, the relevant quantities possess the fohowing values: 



P 



8 (2QoQ^,' - 1) 



5 [bQQliQlf 



6Q0Q0 - 12Q8g'o" - 96Q^(Q^")' - ^QQliQ'^ 



P + Pr 



P + Pt 



a-2Qo 



Ij + 0{x) , 



2 (2QoQo - 1, 
-24Qg(Q^')4 + lOQo(go)' + (^QlQoQo'] + 0{x) 



a-2Qo 



4(2QoQ'o'-l) 



(18i) 

(ISii) 
(ISiii) 



Again it can be noted that the equations of motion become singular where the Ricci 
scalar vanishes. As expected, the degree of the singularity has increased in comparison 
to the n = —1 case. The situation regarding energy conditions here, like for n = — 1, 
does not seem promising in the zero tidal force regime. 

Non-zero tidal force: The following figs. [26] - [ST] display the results of the non-zero 
tidal force scenarios. 





(a) gtt{x) 



(b) p(a 





(c) p + Pr 



(d) p + pt 



Figure 26: R ^ contribution, gu decreasing, throat radius=0.05, and varying xq. For 
negative a -2 the graphs should be flipped around the plane z = 0. 



27 



Figure 27: Anisotropy of throat region: R + a-2R ^ contribution, gtt decreasing, throat 
radius=0.05, a;o = 1, and varying a_2. 




Figure 28: R'^ contribution for xq = 1, gtt decreasing, and varying throat radius. For 
negative a-2 the graphs should be flipped around the plane z = 0. The region around 
the singularity has been omitted. 
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Figure 29: R ^ contribution, gu increasing, throat radius=0.05, and varying xq- For 
negative a_2 the graphs should be flipped around the plane z = 0. 




Figure 30: Anisotropy of throat region: R + a-^R ^ contribution, gu increasing, throat 
radius=0.05, xq = 1, and varying a_2. 
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(a) p{x) (b) p + pr (c) P + Pt 



Figure 31: contribution for xq = 1, g^t increasing, varying throat radius. For 

negative a_2 the graphs should be flipped around the plane z = 0. The region around 
the singularity has been omitted. 



We summarize the results for n = —2 in table 3. 



Tal")lo 3: Suumiary of n -2-^ ^ coutributioiis. 



Parameter 
studied 


zero tidal-force 


gtt{x) concave-down 


gtt{x) concave-up 


xq 


Respecting energy condi- 
tions seems difficult at 
best. 


Energy condition re- 
specting model possible 
near throat. For most 
of positive parameter 
space, small values arc 
found, which may be 
dwarfed by negative 
value of Einstein term). 


Energy condition re- 
specting model not 
found near the throat. 


Throat radius 


Respecting energy condi- 
tions seems diflicult at 
best. 


Seems possible to make 
positive contribution to 
energy conditions for 
very small parameter 
space near the singular 
region for q:_2 < 0. 
(Only on one side of the 
singular region.) 


Does not seem possible 
to respect energy condi- 
tions. 


Minimum 

anisotropy for 


Not studied 


Magnitude of anisotropy 

larger than correspond- 
ing pure Einstein case re- 
gardless of a-2- 


Magnitude of anisotropy 

larger than correspond- 
ing pure Einstein case re- 
gardless of a_2. 
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2.1.6 General n 

Finally, for general n we first quote the throat values in the zero tidal- force case. The 
quantities follows: 

- An^QliQ'^f - SQliQ'if - 6Q0Q0 + 1 " 8n^QliQof + 4ng^QoQo" 

+ 4nQl{Q'^f - SnQliQ'^f - 2nQlQ'^' - WQ^Q'oQT ] + 0{x) , (19i) 

+ nQlQ'^' - 2Ql + 6Qo(Qo)' " 2r?Ql{Q'^ f - 4nQg(Q'^')3 + ^q'^ _ 8Q3(Q'o')' 

- 2nQo(Qo)' - ^QliQ'^ f - QlQ'o - 2nQ[]g'o'go" ] + 0{x) , (19ii) 
p~ + p,=2""i^[l-QoQo'] [l-2Q^Q'^Y-' + 0{x). (19iii) 

These conditions represent the contribution to the energy conditions for a particular 
value of (arbitrary) n. The general energy conditions would then constitute the sum 
of these conditions summed over the n values which contribute to the gravitational 
action. 

For the non-zero tidal force scenarios we present the following results in order to 
show the trends as n increases. The graphs (|32|) -(|35|) are plots of the energy conditions 
exactly at the throat x = for gtt{x) concave-down and gtt{x) concave- up respectively. 
Note that from the analyticity of Q{x) in a neighborhood of the throat, if a quantity 
is positive at the throat then there exists a non-zero neighborhood of the throat where 
this quantity is positive, and hence the particular energy condition can be respected. 
In all the following graphs, only results where the energy conditions are positive are 
shown. (Refer to figure captions for details.) 
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Figure 32: i?" contribution at x = for A = 0.05 and q:„ = 1, gtt decreasing, and 
varying xq. In tlie far-left grapli tlie plots correspond to n = — 3 — > 5 (bottom to top). 
The middle graph corresponds n = —3, —2, —1, 2, 3, 4, 5 (bottom to top). The right 
graph corresponds to parameters n = 1, 2, 3, 4, 5 (bottom to top). Values of n not 
displayed are negative (note though that these would be positive if the sign of a„ was 
negative) . 
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Figure 33: i?" contribution at .t = for xq — 1 and a„ = 1, gn decreasing, and 
varying throat radius. In the far-left graph, the graphs correspond to n = —3 +5 from 
bottom to top. In the middle graph, the plots correspond to n = —3, —2, —1, 3, 4, 5 from 
bottom to top (the others being either negative or too close to zero to show on the scale) . 
In the far-right plot, the curves approaching the vertical- axis correspond to n = 1 — >■ 5 
from bottom to top. The curve that does not approach the vertical-axis corresponds to 
n = -2. 
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(a) p{x) 



(b) p + pr 



(c) P + Pt 



Figure 34: i?" contribution at a; = for A = 0.05 and a„ = 1, gu increasing, and varying 
xo- In the far-left graph the plots correspond to n = — 3 ^ 5 (bottom to top). The middle 
graph corresponds n = 0, —3, —2, —1, 2, 3, 4, 5 from bottom to top. (Note that n = 
should give a contribution of zero, as it corresponds to the cosmological constant only 
scenario. The line at 10~^^ ~ is a numerical artefact.) The right graph corresponds 
to parameters n = 0, 1, 2, 3, 4, 5 (bottom to top, again with the n = plot being a 
numerical artefact.) Values of n not displayed are negative (note though that these would 
be positive if the sign of a„ was negative). 




Figure 35: i?" contribution at x = for xq = 1 and Q!„ = 1, gu increasing, and varying 
throat radius. In the far-left graph, the graphs correspond to n = —3 — )■ +5 from bottom 
to top. In the middle graph, the plots correspond to n = —3, —2, —1, 2, 3, 4, 5 from 
bottom to top (the others being either negative or too close to zero to show on the scale). 
In the far-right plot, the curves intersecting the vertical- axis correspond to n = 1 — )■ 5 
from bottom to top. The curves that do not approach the vertical-axis corresponds to 
n = — 1 (bottom) and n = — 3 (top). 



3 CONCLUDING REMARKS 

We have studied the existence of wormholc throats in modified gravity theories with 
gravitational actions of the form S = Jj^iYl^^n R"'y/9 d'^x. This action includes the 

n 
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n = 1 contribution of Einstein gravity. Of specific interest to our studies are energy 
condition violation and the degree of anisotropy required to support the throat, since 
in Einstein gravity alone the existence of throats implies both violation of energy 
conditions as well as the presence of anisotropy. We studied various terms contributing 
to the action separately, in order to discern which terms aid and which terms hinder 
the energy conditions. It would also be worthwhile to consider all terms together, but 
this would be more computationally involved. We have studied both the zero and 
non-zero tidal force solutions. In general, we find that the parameter space for energy 
condition respecting solutions is much larger for positive n than negative n, and that 
larger n increases the energy conditions. This is due to the fact that f{R) gravity has 
an equivalent description as a scalar-tensor theory, where the "scalar field" component 
may violate energy conditions. The positive n sector also allows one to minimize the 
anisotropy required to support a wormhole throat, whereas the negative n sector tends 
to make the required anisotropy much larger when compared to Einstein gravity. It 
may be then that if gravitation is governed by such a Lagrangian, or one which is 
expandable as such, throats may exist which obey the energy conditions violated by 
Einstein gravity alone and require little or perhaps no anisotropy. In light of this, it 
may be possible that if the true gravitational action contains higher positive powers 
of R, then throats in space-time may be more likely than with no such augmentation 
or with negative powers. Since it is unknown (though perhaps unlikely) whether or 
not the topology of space-time can change, we cannot say whether throats would be 
more likely in the early universe (where higher powers may contribute to infiationary 
expansion) than in the late universe (where inverse powers may contribute to late-time 
acceleration) (see, for example [63] for an analysis of modified gravity accommodating 
early and late-time accelerated expansion). However, it is interesting to speculate on 
this. 
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